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ABSTRACT 

Based on CHANDRA and ASCA observations of nearby starburst galaxies and 
RXTE/ASM, ASCA, and MIR-KVANT/TTM studies of high mass X-ray binary 
(HMXB) populations in the Milky Way and Magellanic Clouds, we propose that the 
number and/or the collective X-ray luminosity of HMXBs can be used to measure 
the star formation rate (SFR) of a galaxy. We show that, within the accuracy of the 
presently available data, a linear relation between HMXB number and star formation 
rate exists. The relation between SFR and collective luminosity of HMXBs is non- linear 
in the low SFR regime, Lx oc SFR~ 1,7 , and becomes linear only for sufficiently high 
star formation rate, SFR > 4.5 M© yr _1 (for M* > 8M©). The non-linear L x - SFR 
dependence in the low SFR limit is not related to non-linear SFR-dependent effects in 
the population of HMXB sources. It is rather caused by the fact, that we measure the 
collective luminosity of a population of discrete sources, which might be dominated by 
the few brightest sources. Although more subtle SFR-dependent effects are likely to 
exist, in the entire range of SFRs the data are broadly consistent with the existence 
of a universal luminosity function of HMXBs which can be roughly described as a 
power law with a differential slope of ~ 1.6, a cutoff at Lx ~ few x 10 40 erg s _1 and 
a normalisation proportional to the star formation rate. 

We apply our results to (spatially unresolved) starburst galaxies observed by 
CHANDRA at redshifts up to z ~ 1.2 in the Hubble Deep Field North and show 
that the calibration of the collective luminosity of HMXBs as a SFR indicator based 
on the local sample agrees well with the SFR estimates obtained for these distant 
galaxies with conventional methods. 
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1 INTRODUCTION 

X-ray observations open a new way to determine the star for- 
mation rate (SFR) in young very distant galaxies. CHAN- 
DRA observations of actively star forming galaxies in our 
vicinity, RXTE/ASM, ASCA, and MIR-KVANT/TTM data 
about high-mass X-ray binary (HMXB) populations in our 
Galaxy and the Magellanic Clouds provide a possibility to 
calibrate the dependence of SFR on the X-ray luminosity 
of a galaxy due to HMXBs. For nearby, spatially resolved 
galaxies for which CHANDRA is able to resolve individual 
X-ray binaries we also have the opportunity to calibrate the 
dependence of SFR on the total number of HMXB sources. 



high mass X-ray binaries, consist of a neutron star (NS) 
or a black hole (BH) accreting from a normal companion 
star. To form a NS or BH the initial mass of the progeni- 
tor star must exceed ~ 8 Mq (Verbunt & van den Heuvel 



1994). The main distinction between LMXBs and HMXBs 



is the mass of the optical companion with a broad, thinly 
populated dividing region between ~ 1 — 5 Mq. This dif- 
ference results in drastically different evolution time-scales 
for low and high mass X-ray binaries and, hence, different 
relations of their number and collective luminosity to the 
instantaneous star formation activity and the stellar con- 
tent of the parent galaxy. In the case of a HMXB, h aving 
a high mass companion, gen erally Af optic( j; > 10 Mq (Ver- 



In ItliP ahgpnrp nf a hrigbt AHW tbp Y-ray pmissmn 



bunt & van den Heuvel 



1994), the characteristic time-scale 



of a galaxy is known to be dominated by the collective 
emissi on of its X-ray binary populations (see e.g. Fabbianc 
(1994)). X-ray binaries, conventionally divided into low and 



is at most the nuclear time-scale of the optical companion 
which does not exceed ~ 2 ■ 10 7 years whereas for a LMXB, 
generally M op u ca i < 1 Mq , it is of the order of • 



10 10 years. 
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HMXBs were first recognised as short-living objects fed 
by the gas supply of a massive star as a result of the dis- 
covery of Cen X-3 as an X-ray pulsar by UHURU, in a bi- 
nary system with an opt ical companion of more than 17 



recently eliminated by Ranalli et al. (2002), who extended 



Mq (Bchreier et al. 1972), and the localisation and mass 



estimation of the Cyg X-l BH due to a soft/hard sta te tran- 



sition occurring simultaneously with a radio flare (Tanan- 



baum et al. 1972), and foll owing optical observations of a 
bright massive counterpart (Bolton 1972| ; Lyutyict al. 1973). 
Dynamics of interacting galaxies, e.g. Antennae, provide an 
additional upper limit on the evolution and existence time- 
scale of HMXBs since the tidal tails and wave patterns in 
which star formation is most vigorous are very short-lived 
phenome na, of the order of a crossing time of inte racting 



galaxies ( Toomre fc Toomre 



1972 



Eneev et al 



1973) 



The prompt evolution of HMXBs makes them a poten- 
tially good t racer of the very rec ent star formation activity 
in a galaxy (Sunyaev et al. 197S) whereas, due to slow evo- 
lution, LMXBs display no direct connection to the present 
value of SFR. LMXBs rather are connected to the total stel- 
lar content of a galaxy determined by the sequence of star 
forma tion episodes experienced by a galaxy d uring its life- 



time flGhosh fc White] |200l| ; 
200^). 



Ptak et al 



2001 



Grimm et al 



Several calibration methods are employed to obtain 
SFRs using UV, FIR and radio flux from distant galaxies. 
Many of these methods rely on a number of assumptions 
about the environment in the galaxy and suffer from vari- 
ous uncertainties, e.g. the influence of dust, escape fraction 
of photons or the shape of the initial mass function (IMF). 
An additional and independent calibrator might therefore 
become a useful method for the determination of SFR. Such 
a method, based on the X-ray emission of a galaxy, might 
circumvent one of the main sources of uncertainty of conven- 
tional SFR indicators - absorption by dust and gas. Indeed, 
galaxies are mostly transparent to X-rays above about 2 
keV, except for the densest parts of the most massive molec- 
ular clouds. 

The existence of various correlations between X-ray and 
optical/far- infrared properties of galaxies has been noted 
and studied in the past. Based on Einstein observations of 
normal and starburst galaxies from the I RAS Bright Galaxy 



Sample, [Griffiths fc Padovani| ( |l990| ) and 



David et al 



(1992) 



found correlations between the soft X-ray luminosity of a 
galaxy and its far-infrared and blue luminosity. Due to the 
limited energy range (0.5-3 keV) of the Einstein observatory 
data one of the main obstacles in quantifying and interpret- 
ing these correlations was proper accounting for absorption 
effects and intrinsic spectra of the galaxies which resulted in 
considerable spread in the derived power law indices of the 
X-ray - FIR correlations, ~ 0.7 — 1.0. Moreover, supernova 
remnants are bright in the soft band of the Einstein obser- 
vatory. CHANDRA, however, is able to distinguish SNRs 
from other sources due to its sensitivity to harder X-rays. 
Although the X-ray data were not sufficient to discriminate 
between contribut ions of different classes of X-ray sources, 
David et al. (1992) suggested that the existence of such cor- 



relations could be understood with a two component model 
for X-ray and far-infrared emission from spiral galaxies, con- 
sisting of old and young populations of the objects having 
different relations to the current star formation activity in 
a galaxy. The uncertainty related to absorption effects was 



these studies to the harder energy band of 2-10 keV based 
on BeppoSAX and ASCA data. In particular, they found a 
linear correlation between total X-ray luminosity of a galaxy 
and both radio and far-infrared luminosities and suggested 
that the X-ray emission might be directly related to the cur- 
rent star formation rate in a galaxy and that such a relation 
might also hold at higher redshifts. 

The main surprise of the study presented here is that 
in the low SFR regime the relation between SFR and col- 
lective luminosity of HMXBs is non-linear, Lx oc SFR~ 17 , 
and becomes linear only for sufficiently high star formation 
rates, when the total number of HMXB sources becomes 
sufficiently large. The non-linear Lx — SFR dependence is 
caused by the fact that we measure the collective luminosity, 
that strongly depends on the brightest sources, of a popula- 
tion of discrete sources. We give a qualitative and approx- 
imate analytic treatment of this (purely statistical) effect 
below and will discus s it in more detail in a separate paper 



ilfanov et al 



2002) 



There are, however, two main obstacles to use the X- 
ray luminosity of a galaxy as a SFR indicator. Firstly, if an 
active nucleus (AGN) is present in a galaxy it can easily 
outshine HMXBs in X-rays. In principle, the presence of an 
AGN component might be identified and, in some cases sep- 
arated, due to different X-ray spectra of an AGN and X-ray 
binaries, provided a sufficiently broad band energy coverage. 
Secondly, there is the dichotomy into LMXBs and HMXBs 
which both have somewhat similar spectra that also could 
probably be distinguished provided sufficiently broad band 
coverage and sufficient signal-to-noise ratio. To estimate the 
SFR one is interested only in the luminosity of HMXBs, 
therefore the LMXB contribution needs to be subtracted. 
This could, in principle, be done based on an estimate of 
the stellar mass of a galaxy. The results of a study of the 



X-ray binary population of our Milky Way (Grimm et al 



2002) and knowledge of the Galactic SFR allow to estimate 
at which point the emission of HMXBs dominates the emis- 
sion of a galaxy. This obviously depends on the ratio of SFR 
to stellar mass of a galaxy. We found, that roughly at a ratio 
of ~ 0.05 Mq yr _1 per 10 10 Mq of total dynamical mass, or 
~ 0.5 Mq yr _1 per 10 10 Mq of stellar mass, the emission of 
HMXB sources begins to dominate the X-ray emission of a 
galaxy (where SFR value refers to a formation rate of stars 
more massive than ~ 5 Mq). It should be emphasised, how- 
ever, that even in the worst case the X-ray luminosity based 
SFR estimate should be able to provide an upper limit on 
the ongoing star formation activity in a galaxy. 

Future observations with present, CHANDRA and 
XMM, and upcoming X-ray missions, Astro-E and especially 
Constellation-X and XEUS, the last having 1 arc sec angular 
resolution and a 100 times larger effective area than CHAN- 
DRA, will permit to get information about the SFR of galax- 
ies from X-rays even at high redshifts. We know from optical 
and radio data that the SFR was m uch higher in galaxies at 
z ~ 2-5 ( |Madau fc Pozzetti| p00o| ). Therefore we could ex- 
pect that in these galaxies the contribution of HMXBs was 
strongly exceeding the contribution of LMXBs. 
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Table 1. The primary sample of local galaxies used to study the luminosity function of HMXB sources. 



Source 


Hubble 


distance'- 1 '-' 


SFR' C ' 


^-dyn arnica I 


Ref.W 


SFR/M 


N(L > 




Lx, total 


Ref.' e ' 




typeM 


[Mpc] 


[M yr" 1 ] 


[1O 1O M ] 




[10- 10 yr -1 ] 


2 ■ 10 38 erg s 




[10 39 erg s" 1 ] 




N3256 


Sb(s) pec 


35.0 


44.0 


5.0 


(i) 


8.8 




12 


128 


(1) 


Antennae 


Sc pec 


19.3 


7.1 


8.0 


(i) 


0.9 




27 


49 


(2) 


M 100 


Sc(s) 


20.4 


4.8 


24.0 


(ii) 


0.2 




5 


10 


(3) 


M 51 


Sbc(s) 


7.5 


3.9 


15.0 


(iii) 


0.26 




15 


16 


(4) 


M 82 


Amorph 


5.7 


3.6 


1.0 


(iv) 


3.6 




12 


23 


(5) 


M 83 


SBc(s) 


3.8 


2.6 


15.4 


(v) 


0.17 




2 


0.14 


(6) 


N4579W 


Sab(s) 


23.5 


2.5 










5 


26 


(7) 


M 74 


Sc(s) 


12.0 


2.2 


14.3 


(vi) 


0.15 




8 


14 


(8) 


Circinus' 9 ' 


Sb 


3.7 


1.5 


2.2 


(v) 


0.73 




6 


5 


(9) 


N4736 


RSab(s) 


4.5 


1.1 


7.0 


(v) 


0.16 




4 


4 


(6) 


N4490 


Scd pec 


8.6 


1.0 


2.3 


(vii) 


0.43 




2 


1.2 


(10) 


N1569 


Sm 


2.1 


0.17 


0.03 


(viii) 


5.6 







0.2 


(11) 


SMC 


Im 


0.06 


0.15 


0.2 


(ix) 


0.75 




1 


0.4 


(12) 


Milky Way 


SAB(rs)bc 


0.25 


54 


(x) 


0.005 





0.2 


(13) 



Mfrom 

(''assuming Hq 



Sandage 



k Tammann 
70 km s -i 



M adopted SFR value ■ 
W References for masses: 



( |l98C|) 

Mpc - 1 and using velocities from Sandage 



Tammann (1980). 



Huchtmcicr fc Richter 



Reakes 




References for X 



Griffiths ct al 



Roberts et al 



• obs ervations: (1 
(6) fcoria fc Wi 
(11) 




fgg§ 



(2) Zeza s 



et al 



_.. ,. d2002[). (3)|Kaarc1 (2001), (4) Icrasliima fc Wilson| (BOOg] ), (5) 

(7) lEraclcous ct al.| (EjrM f 8) poria fc Kont j fflOp . ( 9 ) Smith fc Wilsorj ( |200l| ), (10) 

( |2002| ) 



Martin ct ah] ( |2002j ), (12) |Yokogawa ct at (200C), (13) |Grimm~ct al 



(/J We were not able to obtain a mass value for this source, but according to the rotation curve it is not more massive than the Milky 



Way ( Gonzalez Delgado fc Perez 190o| ) 



Hubble type and velocity taken from ie Vaucouleurs et al. (1991) 



Table 2. The secondary sample of local galaxies used to complement the primary sample in the analysis of the Ljf-SFR relation. 



Source 


Hubble 


distance' 6 ' 


SFR' C ' 


^-dynamical 


Ref.' d ' 


SFR/M 


Lx,total 


Ref.( e ' 




typeM 


[Mpc] 


[M a yr" 1 ] 


[1O 1O M ] 




[lO-Wyr- 1 ] 


[10 39 erg s" 1 ] 




N3690 


Spec 


44.3 


40.0 








220 


(1) 


N7252 


merger 


68.3 


7.7 


4.0 


(i) 


1.9 


94.6 


(2) 


N253 


Sc(s) 


4.2 


4.0 


7.3 


(ii) 


0.55 


5.1 


(3) 


N4945 


Sc 


3.9 


3.1 


9.3 


(ii) 


0.41 


8.9 


(4) 


N3310 


Sbc(r)pec 


15.3 


2.2 


2.0 


(iii) 


1.1 


49.0 


(1) 


N891 


Sb 


11.1 


2.1 


24.0 


(iv) 


0.09 


31.0 


(1) 


N3628 


Sbc 


10.3 


1.6 


16.0 


(v) 


0.1 


13.9 


(5) 


IC342<s) 


Scd 


3.5 


0.48 


11.8 


(ii) 


0.04 


0.9 


(1) 


LMCW 




0.05 


0.25 


0.5 


(vi) 


0.5 


0.34 


(6) 



( a )from Sandage > 
( 6 ) assuming Ho 

(c) 



Tammann 
70 km s -1 



( |l98C|) . 

Mpc - 1 and using velocities from Bandage 



TammannW 198C ) 



adopted SFR value - fro m the last co l umn o f Ta ble | jj or computed from data of [Condon et ah (|l99C|) and Moshir et al. ( [199^) 



erence s for 



X-ra,v ohservatioTis: (TT [Jed a et a,l. X !00j W2) 



Strickland ct al. (2001), (6) Gtrimm ct al 



Re ferenc es for mass es: (i) Li'pari ct a ll (BOO ph, (ii) Huchtmeier & Richter^ (M988h, (iii) 
( 1983j), (v) |5age| (|l993|), (vi) |Feitzingcr| 



Galletta & Recillas-Cruz 



{l982i7Tiv)Pah ca11 



Awaki et al 



( |2002j ), (3) ^.ephaeli fc Grubcr| ( [2002; ), (4) ^church ct al 



'We were not able to obtain a mass value for this source, but the very high SFR ensures for any reasonable mass domination by 
HMXBs. 

Velocity and Hubble type from 



Karachentscv ct al 



CO For LMC data see discussion in the text. 



(1997) 



2 SAMPLE OF GALAXIES 

The list of galaxies used in the following analysis is given in 
Table [I] and Table ^ along with their Hubble type, distances 
and other relevant parameters. 

As our primary sample of local galaxies, used to 
study the HMXB luminosity function and to calibrate 
the L^-SFR relation, we chose a number of nearby late- 



type/starburst galaxies observed by CHANDRA. We based 
our selection primarily on two criteria. Firstly, we selected 
galaxies that can be spatially resolved by CHANDRA suffi- 
ciently well so that the contribution of a central AGN can be 
discriminated and the luminosity functions of the compact 
sources can be constructed without severe confusion effects. 
We should note, however, that for the most distant galaxies 
from our primary sample (e.g. NGC 3256), the probability 
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of source confusion might become non-negligible. Secondly, 
we limited our sample to galaxies known to have high star 
formation rates, so that the population of X-ray binaries 
is dominated by HMXBs and the contribution of lo w m ass 
X-ray binaries can be safely ignored (see subsection 2.5 for 
more detailed discussion). 

In order to probe the HMXB luminosity function in the 
low SFR regime, we used the results of the X-ray binar y 
population study in the Milky Way by Grimm et al. (2002), 



based on RXTE/ASM observations and the luminosity func- 
tion of high mass X-ray b inaries in the Small Magellanic 
Cloud obtained by ASCA ( |Yokogawa et al.|fc000| ). 

The galaxies from our primary sample are listed in Ta- 
ble @. 

In addition, in order to increase the local sample, we 
selected galaxies observed by other X-ray missions, mainly 
ASCA, for which no luminosity function is available, but a 
total flux measurement. The selection was based on the re- 
quirement that no AGN-related activity had been detected 
and the SFR to total mass ratio is sufficiently high to neglect 
the LMXB contribution. These galaxies were used to com- 
plement the primary sample in the analysis of the Lx-SFR 
relation. They are listed in Table |^. 

Finally, in order to study the Lx-SFR relation in dis- 
tant galaxies at redshifts of z ~ 0.2 — 1.3 we used a number of 
galaxies detected by CHANDRA in the Hubble Deep Field 
North, see Table ^. The selection criteria are similar to those 
applied to the local sample and are described in more detail 
in Sec. 3.6. 



2.1 Distances 

To estimate X-ray luminosity and star formation rate, which 
is also based on flux measurements in different spectral 
bands, and compare these values for different galaxies it 
is necessary to have a consistent set of distances. For the 
galaxies from our sample, given in Tables [l] and |^ cosmolog- 
ical effects are not i mportant. The distance s were calculated 
using velocities from Bandage & Tammann ( 1980[ ) corrected 
to the centre of mass of the Local Group and assuming a 
Hubble constant value of Ho = 70 km s~ Mpc - . The dis- 
tances are listed in Table |l| and Table |2| Note that these 
distances might differ from the values used in the original 
publications on the X-ray luminosity functions and SFRs. 



2.2 X-ray luminosity functions 

For the X-ray luminosity functions we used published results 
of Chandra observations of late-type/starburst galaxies. Ref- 
erences to the original publications are given in Table ^ and 
Table |^. The luminosities were rescaled to the distances de- 
scribed in the previous subsection. Note that, due to this 
correction, the total X-ray luminosities and luminosities of 
the brightest sources might differ from the numbers given 
in the original publications. The complete set of luminosity 
functions for all objects from the primary sample (Table fw 
is plotted in Fig. [i] (left panel). 

One of the most serious issues important for the fol- 
lowing analysis is the completeness level of the luminosity 
functions which is obviously different for different galaxies, 
due to different exposure times and distances. In those cases 



when the completeness luminosity was not given in the orig- 
inal publication, we used a conservative estimate based on 
the luminosity at which the luminosity function starts to 
flatten. 

Due to the relatively small field of view of Chandra 
and sufficiently high concentration of X-ray binaries in the 
central parts of the galaxies the contribution of foreground 
and background objects ca n be neglected for the pur pose 
of o ur analysis ( e g. M 51 ( [Terashima fc Wilson||2002[ ), M 
83 ( Soria fc Wu 2002 )). Two of the galaxies in our sam- 
ple - Circinus and NGC 3256 - are located at a Galactic 
latitude of \bn\ < 20°. In these cases the contribution of 
foreground optical stars in the Galaxy that are bright in 
X-rays can be discriminated based on the softness of their 
spectra. Extrapolating t he luminosity funct ion of X-ray bi- 
naries in the Milky Way ( Grimm et al. 2002 ) the probability 
can be estimated of occurence of a foreground source due to 
an unknown Galactic X-ray binary with a flux exceeding the 
sensitivity limit of the corresponding Chandra observations. 
For the Chandra field of view this probability is less than 
~ 10 -3 and therefore can be neglected. 

The luminosities of the compact sources were derived in 
the original publications in slightly different energy bands, 
under different assumptions about spectral shape, and with 
different absorption column densities. Although all these as- 
sumptions affect the luminosity estimates, the resulting un- 
certainty is significantly smaller than those due to distance 
uncertainty and uncertainties in the star formation rate esti- 
mates. Moreover, in many cases, due to insufficient statistics 
of the data an attempt to do corrections for these effects 
could result in additional uncertainties, larger than those 
due to a small difference in e.g. energy bands. Therefore we 
make no attempt to correct for these differences. It should 
be mentioned however, that the most serious effect, up to a 
factor of a few in luminosity might be connected with intrin- 
sic abso rption for the sour ces embedded in dense starforming 
regions ( Zezas et al. 2002 ). Appropriately accounting for this 
requires information about these sources, which is presently 
not available. 

All the luminosity functions with exception of the Milky 
Way are "snapshots" of the duration of several tens of kilo- 
seconds. On the other hand, similar to the Milky Way, com- 
pact sources in other galaxies are known to be variable. E.g. 
NGC 3628 is dominated by a s ingle source, that is kn own to 
vary by about a factor of 30 ( [Strickland et al, 2001). This 
may affect the shape of the individual luminosity functions. 
It should not however affect our conclusions, since in the 
high SFR regime they are based on the average properties 
of sufficiently many galaxies. As for the low SFR regime, the 
Milky Way data are an average of the RXTE/ASM obser- 
vations over four years therefore the contribution of "stan- 
dard" Galactic transient sources is averaged out. 



2.3 Star formation rate estimates 

One of the main uncertainties involved is related to the SFR 
estimates. The conventional SFR indicators rely on a num- 
ber of assumptions regarding the environment in a galaxy, 
such as dust content of the galaxy, or the shape of the ini- 
tial mass function (IMF) . Although comparative analysis of 
different star formation indicators is far beyond the scope of 
this paper, in order to roughly assess the amplitude of the 
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Table 3. The star formation rates for the galaxies from the local sample, measured by different SFR indicators. 



Flux units 
References: (a) 



Source 




Fluxes 








SFRs 








UV 1 


Ha 2 


FIR 3 


radio 4 


Reference 


UV 


Ha 


FIR 


radio 


adopted SFR 


N3256 




0.33 


4.68 




(a) 




5.3 


31.0 












7.1 




00 






47.0 












8.2 




(c) 






54.0 




44 


N4038/9 




1.62 






(d) 




7.9 








(Antennae) 


3.22 


1.36 


4.50 


10.90 


(e) 


9.2 


6.7 


9.0 


9.1 










2.30 




(f) 






4.6 




7.1 


M 100 




0.81 






(d) 




4.5 










3.07 


0.72 


3.36 




(e) 


9.8 


3.9 


7.5 












1.48 




(f) 






3.3 




4.8 


M 51 


15.4 


3.45 


14.7 


8.62 


(e) 


6.6 


2.6 


4.5 


1.1 








4.68 






(d) 




3.5 












2.81 






(g) 




2.1 






3.9 


M 82 




6.17 


52.0 




(h) 




2.6 


9.1 










9.12 






(d) 




3.9 










1.46 


9.98 


112.0 


76.70 


(e) 


0.4 


1.3 


19.6 


5.6 










53.0 




(f) 






9.2 




3.6 


M 83 




13.50 






0) 




2.6 












0.45 






(j) 




0.1 










32.4 


12.20 


34.2 




(e) 


3.6 


2.3 


2.7 




2.6 


N4579 




0.36 






(i) 




2.6 












0.32 






(d) 




2.4 














0.43 




(f) 






1.3 




2.5 


M 74 




1.23 






(d) 




2.3 










6.85 


1.25 


2.92 




(e) 


7.6 


2.4 


2.3 










1.51 






(g) 




2.9 














1.59 




(f) 






1.2 




2.2 


Circinus 






22.3 




(c) 






1.6 










9.5 


16.5 




00 




1.7 


1.2 




1.5 


N4736 




5.37 






(d) 




1.6 
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Figure 1. Left: The luminosity functions of compact X-ray sources in nearby galaxies from the primary sample obtained by CHANDRA 
and listed in Table [l| The luminosity functions are plotted assuming the distances from Table [lj Right: The luminosity functions for the 
same galaxies scaled by the ratio of their star formation rate to the SFR of Antennae. The luminosity functions in the right panel are 
plotted only above their corresponding completeness limits. It is clear that despite large differences in the SFRs (by a factor of ~ 40 — 50) 
the scaled luminosity functions occupy only a narrow band in the 7V(> L) — L plane. 



uncertainties in the SFR estimates we compared results of 
different star formation indicators for each galaxy from our 
sample with special attention given to the galaxies from the 
primary sample. For all galaxies from the primary sample we 
found at least 3 different measurements of star formation in- 
dicators in the literature, namely UV, Ha, FIR or thermal 
radio emission flux. The data along with the corresponding 
references are listed in Table ^. 

In order to convert the flux measurements to star forma- 
tion rates we use the result of an empirical crosscalibration 



of sta r formation rate indicators by Rosa-Gonzalez et al 
(2002^. The calibration is based on the canonical formulae 



by [Kcnnicutt ( 1998 ) and takes into account dust /extinction 
effects. We used the following flux-SFR relations: 



SFR Ha = 1.1-10" 
SFRuv = 6.4 ■ 1CT 28 
SFR FIR = 4.5 • 10" 
SFR radlo = 1.82 ■ 10" 



■ L Ha [erg s~ ] 
• Luv [erg s"" 1 Hz -1 ] 
44 ■ L FIR [eig s" 1 ] 
■ v dh z ■ L v [erg s' 1 Hz^ 1 ] 



(1) 
(2) 
(3) 
(4) 



The last relation is from 



Condon 



(1992) and applies only 
to the thermal radio emission, originating, presumably, in 
hot gas in HII regions associated with star formation (as we 
used thermal 1.4 GHz flux estimates from Bell & Kennicutt 



( pooip ). 

The above relations refer to the SFR for stars more mas- 
sive than ~ 5 Mq. The total star formation rate, including 
low mass stars, could theoretically be obtained by extrap- 
olating these numbers assuming an initial mass function. 
Obviously, such a correction would rely on the assumption 
that the IMF does not depend on the initial conditions in 
a galaxy and would involve a significant additional uncer- 



tainty. On the other hand, this correction is not needed for 
our study as the binary X-ray sources harbour a compact 
object - a NS or a BH - which according to the modern 
picture of stellar evolution can evolve only from stars with 
initial masses exceeding ~ 8 Mq. The SFR correction from 
M > 5 M to M > 8 M is relatively small (~ 20 per 
cent) and, most importantly, due to the similarity of the 
IMFs for large masses it is significantly less subject to the 
uncertainty due to poor knowledge of the slope of the IMF. 
Thus, for the purpose of our study it is entirely sufficient to 
use the relations (|l|)-(§) without an additional correction. 
In the following, the term SFR refers to the star formation 
rate of stars more massive than ~ 5 Mq. 



Since the relations (|l|)-(|4|) are based on the average 
properties of star forming galaxies there is considerable scat- 
ter in the SFR estimates of a galaxy obtained using different 
indicators (Table ^). On the other hand, the SFR estimates 
based on different measurements of the same indicator are 
generally in a good agreement with each other. A detailed 
study, which SFR indicator is most appropriate for a given 
galaxy is beyond the scope of this paper. Therefore, we re- 
lied on the fact that for all galaxies from our primary sample 
there are more than 3 measurements for different indicators. 
For each galaxy we disregarded the estimates significantly 
deviating from the majority of other indicators, and aver- 
aged the latter. The final values of the star formation rates 
we used in the following analysis are summarised in the last 
column of Table |5| 
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Figure 2. Contributions of LMXBs and HMXBs to the observed 
luminosity function for NGC 4736 (thin solid histogram), hav- 
ing smallest SFR to total mass ratio in the primary sample. The 
upper thick grey histogram corresponds to the contribution of 
HMXBs scaled from the Milky Way HMXB luminosity function 
by the ratio of the SFRs. The lower dotted histogram is the Galac- 
tic LMXB luminosity function scaled by the ratio of the total 
masses. Total masses and SFRs are given in Table [lj. 

2.4 Contribution of a central AGN 

As mentioned in Sec. ^ the emission of a central AGN can 
easily outshine the contribution of X-ray binaries. However, 
due to the excellent angular resolution of CHANDRA it is 
possible to exclude any contribution from the central AGN 
in nearby galaxies. In our primary sample a central AGN 
is present in the Circinus galaxy and NGC 4579 for which 
the point source associated with the nucleus of the galaxy 
was excluded from the luminosity function. Also NGC 4945 
is a case where there is contribution to the X-ray emission 
from an AGN. However the AGN is heavily obscured and the 
emission be low about 10 keV of the AGN negligible ( schurcb 
et aDpOoi). 



2.5 Contribution of LMXBs 

Due to the absence of optical identifications of a donor star 
in the X-ray binaries detected by CHANDRA in other galax- 
ies, except for LMC and SMC, there is no obvious way to 
discriminate the contribution of low mass X-ray binaries. On 
the other hand the relative contribution of LMXB sources 
can be estimated and, as it was mentioned above, it was one 
of the requirements to minimise the LMXB contribution, 
that determined our selection of the late-type/starburst 
galaxies. 

Due to the long evolution time-scale of LMXBs we ex- 
pect the population of LMXB sources to be roughly propor- 
tional to the stellar mass of a galaxy, whereas the population 
of short-living HMXBs should be defined by the very recent 
value of the star formation rate. Therefore the relative im- 



portance of LMXB sources should be roughly characterised 
by (inversely proportional to) the ratio of star formation rate 
to stellar mass of a galaxy. Since the determination of stellar 
mass, especially for a starburst galaxy, is very difficult and 
uncertain we used values for the total mass of a galaxy esti- 
mated from dynamical methods and assumed that the stel- 
lar mass is roughly proportional to the total mass. To check 
our assumption we compare the dynamical mass with the 
K band luminosity for galaxies for which, first, enough data 
exist to construct a growth curve in the K band and, second, 
for which an extrapolation to t he total K band flux c an be 
made following the approach of Spinoglio et al. ( |lOU5| ). The 
number of galaxies is small, the sample consists of M 74, M 
83, NGC 4736 and NGC 891, and the uncertainties associ- 
ated with this approach are big, i.e. of order a factor 3. But 
within this uncertainty there is a correlation between the K 
band luminosity and the dynamical mass estimate. However, 
due to the more abundant data for and higher accuracy of 
dynamical masses we do not use stellar mass estimates based 
on K band luminosities in the following. The values of the 
total dynamical mass, corresponding references, and the ra- 
tios of SFR to total mass are given in Table tl and Table 

I 

The SFR to total mass ratios for late-type galaxies 
should be compared with that for the Milky Way, for which 
the population of s ufficiently luminou s X-ray binaries is 
studied rather well (Glrimm et al. 2002). We know that the 



Milky Way, having a ratio SFR/M d »n ~ 5 • 10" 13 yr" 1 , 
or SFR/M ste iiar ~ 5 • 10" 12 yr" 1 , is dominated by LMXB 
sources, HMXB sources contributing ~ 10 per cent to the 
total X-ray luminosity and ~ 15 per cent to the total num- 
ber (above ~ 10 37 erg s" 1 ) of X-ray binaries. As can be seen 
from Table |l], concerning the galaxies for which luminos- 
ity functions are available the minimal value of SFR/M<j yn 
~ 1.5- 10" 11 yr" 1 is achieved for M 74 and NGC 4736, which 
exceeds by a factor of ~ 30 that of the Milky Way. There- 
fore, even in the least favourable case of these two galaxies, 
we expect the HMXB sources to exceed LMXBs by a factor 
of ~ 3 at least, both in number and in luminosity. A more 
detailed comparison is shown in Fig. |^, where we plot the 
expected contributions of LMXBs and HMXBs to the ob- 
served luminosity function for NGC 4736. The luminosity 
function of HMXBs was obtained by scaling the Milky Way 
HMXB luminosity function by the ratio of SFRs of NGC 
4736 to the Milky Way. The LMXB contribution was simi- 
larly estimated by scaling the Milky Way LMXB luminosity 
function by the ratio of the corresponding total masses. As 
can be seen from Fig. the contribution of LMXB sources 
does not exceed ~ 30 per cent at the lower luminosity end 
of the luminosity function. If the fractions of NSs and BHs 
in low mass systems in late-type/starburst galaxies are sim- 
ilar to that in the Milky Way, the contribution of LMXBs 
should be negligible at luminosities above ~ 10 38 erg s , 
corresponding to the Eddington limit of a neutron star, to 
which range most of the following analysis will be restricted. 

For all galaxies from Tables H and ^ the lowest values for 
SFR/M are 4 ■ 10" 12 and 9 ■ 10" H for IC 342 and NGC 891, 
respectively. This means that the contribution of LMXBs 
could make up a sizeable portion of their X-ray luminosity, 
~50 per cent for IC 342 and -25 per cent for NGC 891. 
Therefore their data points should be considered as upper 
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SFR [M /yr] 

Figure 3. Number of sources with a 2-10 keV luminosity exceed- 
ing 2 TO 38 erg s _1 versus SFR for galaxies from Table [lj The figure 
shows a clear correlation between the number of sources and the 
SFR. The straight line is the best fit from a Maximum-Likelihood 
fit, Eq. |H| The vertical error bars were calculated assuming a Pois- 
sonian distribution, the SFR uncertainty was assumed to be 30 
per cent. For M 74 and M 100, whose completeness limit exceeds 
2 • 10 38 erg s _1 the contribution of sources above 2 ■ 10 38 erg 
s" 1 and below the completeness limit was estimated from the 
"universal" luminosity function, Eq. j^. 



limits on the integrated luminosity of HMXBs (shown in 
Fig. as arrows) . 



3 HIGH MASS X-RAY BINARIES AS A STAR 
FORMATION INDICATOR 

As already mentioned, the simplest assumption about the 
connection of HMXBs and SFR would be that the number 
of X-ray sources with a high mass companion is directly pro- 
portional to the star formation rate in a galaxy. In Fig. [j] 
(right panel) we show the luminosity functions of the galax- 
ies from our primary sample scaled to the star formation 
rate of the Antennae galaxies. Each luminosity function is 
plotted above its corresponding completeness limit. It is ob- 
vious that after rescaling the luminosity functions occupy a 
rather narrow band in the log(N)-log(L) plane and seem to 
be consistent with each other within a factor of ~ 2 whereas 
the star formation rates differ by a factor of ~ 40 — 50. 
This similarity indicates that the number/luminosity func- 
tion of HMXB sources might indeed be proportional to the 
star formation rate. This conclusion is further supported by 
Fig. ^ which shows the number of sources with a luminosity 
above 2 • 10 38 erg s -1 versus the SFR. The threshold lu- 
minosity was chosen at 2 ■ 10 38 erg s" 1 to have a sufficient 
number of galaxies with a completeness limit equal or lower 
than that value and, on the other hand, to have a sufficient 
number of sources for e ach individual galaxy. In addition, 
as was discussed in Sec. E.5l this choice of the threshold lu- 
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Figure 4. Comparison of the combined luminosity function of 
M 82, NGC 4579, NGC 4736 and Circinus, having SFRs in the 
range 1-3.5 Mq yr -1 with the Antennae luminosity function (7.1 
Mq yr — 1 ). A Kolmogorov-Smirnov test gives a probability of 15 
per cent that the two luminosity functions are derived from the 
same distribution. See discussion in the text regarding the effect 
of the errors in the distance measurements on the shape of the 
combined luminosity function. 



minosity might help to minimise the contribution of LMXB 
sources. The errors for the number of sources were computed 
assuming a Poissonian distribution. For the SFR values we 
assumed a 30 per cent uncertainty. Although the errors are 
rather big, the correlation of the number of sources with SFR 
is obvious. The slope of the correlation, determined from a 
least-squares fit in the form N oc SFR a , is a — 1.06 ± 0.07, 
i.e. it is consistent with unity. A fit of this correlation with 
a straight line N oc SFR (shown in the figure by solid line) 



N(L > 2 • 10 38 erg s" 1 ) = (2.9 ± 0.23) ■ SFR[M Q yr" 



(5) 



According to this fit we should expect less than 1 source in 
the Milky Way, having a SFR of 0.25 M Q yr -1 , which is in 
agreement w ith the fact that no source above this luminosity 
is observed (Grimm et al. 2002). 



3.1 Universal HMXB Luminosity Function ? 

In order to check the assumption that all the individual lu- 
minosity functions have identical or similar shape with the 
normalisation being proportional to the SFR, we compare 
the luminosity function of the Antennae galaxies, having a 
high star formation rate (~ 7 Mq yr -1 ), with the collective 
luminosity function of galaxies with medium SFRs (in the 
range of ~ 1.0-3.5 M yr" 1 ). For the later we summed the 
luminosity functions of M 82, NGC 4579, NGC 4736 and 
Circinus, having a combined SFR of ~ 8.8 Mq yr -1 . The 
two luminosity functions (shown in Fig. |^) agree very well at 
Lx < 10 39 erg s _1 with possible differences at higher lumi- 
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nosities. In a strict statistical sense, a Kolmogorov-Smirnov 
test gives a 15 per cent probability that the luminosity func- 
tions are derived from the same distribution, thus, neither 
confirming convincingly, nor rejecting the null hypothesis. 
However, it should be emphasised, that whereas the shape 
of a single slope power law luminosity function is not affected 
at all by the uncertainty in the distance, more complicated 
forms of a luminosity function, e.g. a power law with cut- 
off, would be sensitive to errors in the distance determina- 
tion. The effect might be even stronger for the combined 
luminosity functions of several galaxies, located at different 
distances and each having different errors in the distance 
estimate. In the case of a power law with high luminosity 
cut-off, the effect would be strongest at the high luminos- 
ity end and will effectively dilute the cut-off, as probably 
is observed. Therefore, we can presently not draw a defini- 
tive conclusion about the existence of a universal luminosity 
function of HMXBs, from which all luminosity functions of 
the individual galaxies are strictly derived. For instance, sub- 
tle effects similar to the effect of flattening of the luminosity 
function with increase of SFR suggested by [Kilgard et al" 



( t200l| ); |Ghosh fe White] ( |200l| ) 



Ptak et al 



( 2001) can not be 



excluded based on the presently available sample of galaxies 
and sensitivities achieved. We can conclude, however, that 
there is no evidence for strong non-linear dependences of the 
luminosity function on the SFR. 

As the next step we compare the luminosity functions 
of actively star forming galaxies with that of low SFR galax- 
ies. Unfortunately, the X-ray binary population of low SFR 
galaxies is usually dominated by LMXB systems. One of the 
cases in which the luminosity function of HMXB sources can 
be reliably obtained is the Milky Way galaxy, for which all 
sufficiently bright X-ray binaries are optically identified. An- 
other case is the Small Magellanic Cloud, which has a SFR 
value similar to our Galaxy, but is ~ 300 — 500 less massive 
and, correspondingly, has very few, if any, LMXB sources 
(fokogawa et al. 200C). Moreover, the SMC is close enough 
to have optical identifications of HMXBs which makes a dis- 
tinction like in the Milky Way possible. In order to do the 
comparison, we combined the luminosity functions of all ac- 
tively star forming galaxies from our sample with a com- 
pleteness limit lower than 2 • 10 38 erg s _1 - M 82, Antennae, 
NGC 4579, NGC 4736 and Circinus. These galaxies have a 
total SFR of ~ 16 M yr -1 , which exceeds the Milky Way 
SFR (~ 0.25 M yr" 1 ) by a factor of ~ 65. Fig. | shows the 
combined luminosity function of the above mentioned star 
forming galaxies and the luminosity functions of Galactic 
and SMC HMXBs scaled according to the ratios of SFRs. 
Shown in Fig. | by a solid line is the fit to the luminosity 
function of the high SFR galaxies only (see below) , extrapo- 
lated to lower luminosities. It is obvious that the luminosity 
functions of Galactic and SMC HMXBs agree surprisingly 
well with an extrapolation of the combined luminosity func- 
tion of the starburst galaxies. 

Thus we demonstrated that the presently available 
data are consistent with the assumption that the approxi- 
mate shape and normalisation of the luminosity function for 
HMXBs in a galaxy with a known star formation rate can be 
derived from a "universal" luminosity function whose shape 
is fixed and whose normalisation is proportional to star for- 
mation rate. Due to a number of uncertainties involved, the 
accuracy of this approximation is difficult to assess. Based 



on our sample of galaxies we can conclude that it might be 
accurate within ~ 50 per cent. 

In order to obtain the universal luminosity function 
of HMXBs we fit the combined luminosity function of M 
82, Antennae, NGC 4579, NGC 4736 and Circinus using a 
Maximum-Likelihood method with a power law with a cut- 
off at L c — 2.1 ■ 10 40 erg s" 1 and normalise the result to 
the combined SFR of the galaxies. The best fit luminosity 
function (solid line in Figj5J) in the differential form is given 
by: 



dN 



(3.3t£ 



1-1 ^ ■ SFR ■ LJg 1 ' 61 * 0,12 forL<L c , (6) 



where L38 = Z//10 38 erg s" 1 and SFR is measured in units 
of Mq yr -1 . The errors are la estimates for one parameter 
of interest. The rather large errors for normalisation are due 
to the correlation between slope and normalisation of the 
luminosity function, with a higher value of normalisation 
corresponding to a steeper slope. The cumulative form of 
the luminosity function, corresponding to the best values of 
the slope and normalisation is: 



N(> L) = 5.4 • SFR ■ (LI 



210~ 



) 5 



(7) 



According to a Kolmogorov-Smirnov test the data are con- 
sistent with the best fit model at a confidence level of 90 per 
cent. 

As an additional test we checked all individual lumi- 
nosity functions against our best fit using a Kolmogorov- 
Smirnov test. Taking into account the respective complete- 
ness limits, the shapes of all individual luminosity functions 
are compatible with the assumption of a common 'origin'. 
In Fig. ^| we show the individual luminosity functions along 
with the universal luminosity function given by Eq.(^|) with 
the normalisation determined according to the correspond- 
ing star formation rates derived from the conventional SFR 
indicators (Table Q). 

Finally, we construct the differential luminosity func- 
tion combining the data for all galaxies from the primary 
sample, except for NGC 3256 (having somewhat uncertain 
completeness limit) . To do so we bin all the sources above the 
corresponding completeness limits in logarithmically spaced 
bins and normalise the result by the combined SFR of all 
galaxies contributing to a given bin. Such a method has the 
advantage of using all the available data. A disadvantage is 
that due to significantly different luminosity ranges of the 
individual luminosity functions (especially SMC and Milky 
Way on one side and star forming galaxies on the other) un- 
certainties in the conventional SFR estimates may lead to 
the appearance of artificial features in the combined lumi- 
nosity function. With that in mind, we plot the differential 
luminosity function in the right panel of Fig. |E] along with 
the best fit power law from Eq.(^|). 

In order to investigate the influence of systematic un- 
certainties in SFR and distance we performed a Monte-Carlo 
simulation taking into account these two effects. The grey 
area in the right panel of Fig. [| shows the 90 per cent confi- 
dence interval obtained from this simulation. In the simula- 
tion we randomly varied the distances of galaxies, assuming 
the errors on the distance to be distributed according to a 
Gaussian with a mean of and a width of 20 per cent of the 
distance of a galaxy which corresponds to an uncertainty 
in luminosity of ~ 40 per cent. Correspondingly the SFR, 
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Figure 5. Left: Combined luminosity function of compact X-ray sources in the starburst galaxies M82, NGC 4038/9, NGC 4579, NGC 
4736 and Circinus with a total SFR of 16 Mq yr — 1 (histogram above 2 • 10 38 erg s —1 ) and the luminosity functions of NGC 1569 and 
HMXBs in the Milky Way and Small Magellanic Clouds (three histograms below 2 • 10 38 erg s — 1 ). The thin solid line is the best fit 
to the combined luminosity function of the starburst galaxies only, given by Eq. |?j Right: Differential luminosity function obtained by 
combining the data for all galaxies from the primary sample, except for NGC 3256 (see text). The straight line is the best fit to the 
luminosity function of star forming galaxies given by Eq. |i| - the same as in the left hand panel. Note, that due to different construction 
algorithms, the luminosity functions shown in the left and right panels are based on different but overlapping samples of galaxies (see 
discussion in the text). The grey area is the 90 per cent confidence level interval we obtained from a Monte-Carlo simulation taking into 
account uncertainties in the SFR and distances. For details see discussion in the text. 



affected in the same way as the X-ray luminosity by uncer- 
tainties in the distance, was changed. Additionally the SFR 
was randomly varied also assuming a Gaussian error distri- 
bution with a mean of and a width of 30 per cent of the 
SFR, as assumed for Fig. ^[ For the Milky Way we varied in 
each Monte-Carlo run the distance to each HMXB indepen- 
dently with a Gaussian with a mean of and a width of 20 
per cent of the distance. 

Noteworthy is the fact that the luminosity function is 
sufficiently close to a single slope power law in a broad lu- 
minosity range covering more than five orders of magnitude. 
If the absence of significant features is confirmed this allows 
to constrain the relative abundance of NS and BH binaries 
and/or the properties of accreting compact obj ects at su- 
percritical accretion rates (see discussion in Sec. |4.l[ ). 

However, it should be emphasised that there is hardly 
any overlap in the luminosity functions for low and high 
SFR galaxies, as is obvious from Fig. [j] and ^. It happens 
that this gap is around the Eddington luminosity of a NS, 
L/Edd.NS, which should be a dividing line between NS and 
BH binaries. From simple assumptions it would be expected 
that the luminosity functions below LEdd,NS are dominated 
by NS whereas above Lsdd.NS BH binaries should domi- 
nate. This would imply a break in the luminosity function 
around Lsdd,NS because of different abundances of NSs and 
BHs. Due to the uncertainties in SFR measurements it is 
possible that a break, that would theoretically be expected 
around L E dd,NS, could be hidden by this gap. Even upper 
limits (not more than twice) are of importance and could 



give some additional information about the relative strength 
of the two populations of accreting binaries (see discussion 
in Sec. ^). Observations of star forming galaxies with suf- 
ficient sensitivity, i.e. with a completeness limit well below 
10 38 erg s _1 will be able to resolve this question. 

3.2 High Luminosity cut-off 

The combined luminosity function shown in the left panel of 
Fig. [| indicates a possible presence of a cut-off at L c ~ 2T0 40 
erg s _1 . From a statistical point of view, when analysing 
the combined luminosity function of the high SFR galaxies 
only, the significance of the cut-off is not very high, with a 
single slope power law with slope 0.74 for the cumulative 
luminosity function also giving an acceptable fit, although 
with a somewhat lower probability of 54 per cent according 
to a Kolmogorov-Smirnov test. However, an independent 
strong evidence for the existence of a cut-off around ~ few • 
10 40 erg s _1 is provided by the Ljf-SFR relation as discussed 
in the next subsections. 

The existence of such a cut-off, if it is real and if it is a 
universal feature of the HMXB luminosity function, can have 
significant implications to our understanding of the so-called 
ultra-luminous X-ray sources. Assuming that these very lu- 
minous objects are intermediate mass BHs accreting at the 
Eddington limit, the value of the cut-off gives an upper limit 
on the mass of the black hole, ~ 100 Mq . These apparently 
super-Eddington luminosities can also be the result of other 
effects, like a strong magnetic field in NSs which may allow 
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Figure 6. Comparison of the "universal" luminosity function defined by Eq. ^ (thin solid lines) with individual luminosity functions of 
compact X-ray sources in the galaxies from Table hi (histograms). The normalisation of the "universal" luminosity function in each panel 
was calculated using corresponding SFR values from Table [l]. For M83 the luminosity function of the compact sources in the nuclear 
region only is plotted, whereas the normalisation of the "universal" luminosity function was computed using the overall SFR for the 
galaxy. Therefore the thin line should be considered as an upper limit. The dotted lines are fits to the normalisation of the observed 
luminosity functions in the cases where completeness or coverage do not represent the same area as the SFR measurements. 
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radiation to escape without inter acting with the accreting 
material ( Basko fc Surry acv |1976D, emission from a super- 



critical accretion disk ( Shakura fc Sunyacv |l973t |Paczyn , 
sky & Wiita 1980), beamed emission (|King et*l=iL[ 2001 ), or 



the luminosity of t he next brightest sou rce is by a factor of 
~ 30 - 50 smaller flSunyaev etaT| |l99d|). 



the emission of a jet as suggested by Kording et al. (2002) 
Moreover, in BHs in high state radiation is coming from the 
quasi-flat accretion disk where electron scattering gives the 
main contribution to the opacity. It is easy to show that 
the radiation flux perpendicular to the p lane of the disk 



exce eds th e average value by up to 3 times ( Shakura fc Sun 



yaev 



It 73). Also the Eddington luminosity is dependent on 



chemical abundance which allows a twice higher luminosity 
for accretion of helium. These last two effects alone can pro- 
vide a factor of 6 above the canonical Eddington luminosity. 

It should be mentioned that, based on the combined 
luminosity function only, we can not exclude the possibility 
that the cut-off is primarily due to the Antennae galaxies 
which contributes about half of the sources above 10 39 erg 
s _1 and shows a prominent cut-off in its luminosity function. 
On the other hand, further indication for a cut-off is pro- 
vided by the luminosity function of NGC 3256. Conventional 
star formation indicators give a value of SFR of ~ 45 Mq 
yr _1 , however its luminosity function also shows a cut-off 
at ~ 10 40 erg s _1 . Unfortunately, due to the large distance 
(35 Mpc) and a comparatively short exposure time of the 
CHANDRA observation, ~ 28 ks, the luminosity function 
of NGC 3256 becomes incomplete at luminosities shortly 
below the brightest source and therefore does not allow for 
a detailed investigation. 



3.3 Total X-ray luminosity as SFR indicator 

CHANDRA and future X-ray missions with angular reso- 
lution of the order of ~ 1" would be able to spatially re- 
solve X-ray binaries only in nearby galaxies (d< ~ 50— 100 
Mpc). For more distant galaxies only the total luminosity of 
a galaxy due to HMXBs can be used for X-ray diagnostics 
of star formation. 

Fig. in shows the total luminosity of X-ray binaries 
(above 10 erg s _1 ) plotted versus SFR. The galaxies from 
the primary sample (listed in Table |l]) are shown by filled 
circles. The galaxies, for which only total luminosity is avail- 
able (Table W are shown as filled triangles. The luminosities 
of the galaxies from the primary sample were calculated by 
summing the luminosities of individual sources down to the 
completeness limit of the corresponding luminosity function. 
The contribution of the sources below the completeness limit 
was approximately accounted for by integrating a power law 
distribution with slope 1.6 and normalisation obtained from 
the fit to the observed luminosity function. Note, that due 
to the shallow slope of the luminosity function the total lu- 
minosity depends only weakly on the lower integration limit. 

As an additional data point we take luminosity and SFR 
for the Large Ma gellanic Cloud. The S FR is similar to the 
Milky Way SFR ( |Holtzman et al.| [l999| ). Since no luminos- 
ity function is presently available for LMC we estimated its 
integrated X-ray luminosity as a sum of the time averaged 
luminosities of the three brightest HM XB sources (LMC X - 



2002) 



1, X-3, X-4) as measured by ASM (Grimm et al 
£2-10 kcv ~ 3.4 ■ 10 38 erg s _1 . Contribution of the weaker 
sources should not change this estimate significantly, since 



3.4 Theoretical L X -SFR relation 

At first glance, the relation between collective luminosity 
of HMXBs and SFR can be easily derived integrating Eq. 
([]) for the SFR dependent luminosity function. Therefore, 
as the population of HMXB sources in a galaxy is directly 
proportional to SFR, one might expect that the X-ray lu- 
minosity of galaxies due to HMXB, Lx, should be linearly 
proportional to SFR. However this problem contains some 
subtleties related to the statistical properties of the power 
law luminosity distribution of discrete sources which appear 
not to have been recognised previously (at least in astro- 
physical context). The difference between the most prob- 
able value of the total luminosity of HMXB sources in a 
galaxy (the mode of the distribution) and the ensemble av- 
erage value (expectation mean, obtained by integrating Eq. 
(^)) results in the non-linear Lx-SFR dependence in the 
low SFR regime. As this effect might be of broader general 
interest and might work in many different situations related 
to computing/measuring integrated luminosity of a limited 
number of discrete objects, we give it a more detailed an d 
rigorous discussion in a separate paper (Gilfanov et al, 2002), 



and restrict the discussion here to only a brief explanation 
and an approximate analytical tr eatment. A somewhat s im- 
ilar problem was considered by Kalogera et al. ( 2001 ) in 
the context of pulsar counts and the faint end of the pulsar 
luminosity function. 

For illustration only, let us consider a population of dis- 
crete sources with a Gaussian luminosity function. As is well 
known, in this case the sum of their luminosities - the inte- 
grated luminosity of the parent galaxy, also obeys a Gaus- 
sian distribution for which the mean luminosity and disper- 
sion can be computed straightforwardly. An essential prop- 
erty of this simple case is that for an ensemble of galaxies, 
each having a population of such sources, the most probable 
value of the integrated luminosity of an arbitrarily chosen 
galaxy (the mode of the distribution) equals to the mean lu- 
minosity (averaged over the ensemble of galaxies) . The situa- 
tion might be different in the case of a population of discrete 
sources with a power law (or similarly skewed) luminosity 
function. In this case an ensemble of galaxies would have a 
non-Gaussian probability distribution of the integrated lu- 
minosity. Due to skewness of the probability distribution in 
this case, the most probable value of the integrated lumi- 
nosity of an arbitrarily chosen galaxy does not necessarily 
coincide with the mean value (the ensemble average). The 
effect is caused by the fact that depending on the slope of 
the luminosity function and its normalisation the integrated 
luminosity of the galaxy might be defined by a small number 
of brightest sources even when the total number of sources 
is large. Of course, in the limit of large number of sources 
in the high luminosity end of the luminosity function the 
distribution becomes asymptotically close to Gaussian and, 
correspondingly, the difference between the most probable 
value and the ensemble average vanishes. In this limit the 
relation between the integrated luminosity of HMXBs and 
SFR can be derived straightforwardly integrating Eq.(^|) for 
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L c = 2.1 ■ 10 40 erg s _1 : 



.7 • 10 39 ■ SFR[M e yr _1 ] erg s" 1 



(8) 



It should be emphasised that the ensemble average inte- 
grated luminosity (i.e. averaged over many galaxies with 
similar SFR) is always described by the above equation, in- 
dependent of the number of sources and shape of the lu- 
minosity function. This equality is maintained due to the 
outlier galaxies, whose luminosity exceeds significantly both 
the most probable and average values. These outlier galax- 
ies will result in enhanced and asymmetric dispersion in the 
low SFR-regime. 

The following simple consideration leads to an approx- 
imate analytical expression for the most probable value of 
the integrated luminosity. Assuming a power law luminosity 
function dN/dL = A ■ SFR ■ L~ a with 1 < a < 2, one 
might expect, that the brightest source would most likely 
have a luminosity L ma x close to the value ~ L\ such that 
N(> Li) ~ 1, i.e. 



Li oc SFR— 



(9) 



In the presence of a cut-off L c in the luminosity function, 
the luminosity of the brightest source, of course, can not 
exceed the cut-off luminosity: L max = min(Li,L c ). The 
most probable value of the total luminosity can be com- 
puted integrating the luminosity function from L m in to 
Lmax = min(Li, L c ): 

rmin(L 1 ,L c ) jju 

' " (10) 



-L dL, 



which leads to 

A ■ SFR 



J total 



dL 



min(Li, L c ) 2 



(11) 



for 1 < a < 2 and Li,L c » L m i n . 

Obviously there are two limiting cases of the Lx-SFR 

dependence of the total luminosity on SFR, depending on 

the relation between L c and Li, i.e. on the expected number 

of sources in the high end of the luminosity function, near 

its cut-off. In the limit of low SFR (small number of sources) 

Li < L c and the luminosity of the brightest source would 

i 

increase with SFR: L mal ~ L\ oc SFR"- 1 . Therefore the 
Lx-SFR dependence might be strongly non-linear: 



Ltotai oc SFR 7 



(12) 



e.g. for a — 1.5 the relation is quadratic Ltotai oc SFR . For 
sufficiently large values of SFR L\ > L c , i.e. N(> L c ) > 1 
implying a large number of sources in the high luminosity 
end of the luminosity function and, correspondingly, Gaus- 
sian probability distribution of the integrated luminosity. In 
this case L max ~ L c — const and does not depend on SFR 
anymore and the dependence is linear, in accord with Eq.(|^). 

Importantly, the entire existence of the linear regime 
in the Lx-SFR relation is a direct consequence of the exis- 
tence of a cut-off in the luminosity function. For a sufficiently 
flat luminosity function, 1 < a < 2, the collective luminos- 
ity of the sources grows faster than linear because brighter 
and brighter sources define the total luminosity as the star 
formation rate increases. Only in the presence of the max- 
imum possible luminosity of the sources, L c (for instance 
Eddington limit for NSs) the regime can be reached, when 
7V(> L c ) becomes larger than unity and subsequent increase 
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Figure 7. The Lx -SFR relation. The filled circles and triangles 
are nearby galaxies from Table [l] (primary sample) and Table 
^ (secondary sample), the open circles are distant star forming 
galaxies from the HDF North and Lynx field. The arrows are the 
upper limits for the X-ray luminosity due to HMXBs for IC 342 
and NGC 891. The thick solid line shows the expected relation 
between SFR and the most probable value of the total luminosity 
computed for the best fit parame ters of the HMXB lu minosity 
function (exact calculation, from Gilfanov et al. (2002)). Note, 



that in the low SFR regime the probability to find a galaxy below 
the solid curve is ~ 10—15 per cent. The shaded area shows the 68 
per cent confidence region including both intrinsic variance of the 
Lx— SFR relation and uncertainty of the best fit parameters of 
the HMXB luminosity function (Eq. (^)). The dashed line shows 
the linear Lx-SFR relation given by Eq. (H). 



of the star formation rate results in the linear growth of the 
total luminosity. The latter, linear, re gime of the Lx-SFR 

~ (pool) 



Ranalli et al 



relation was studied independently by 
based on ASCA and Beppo-SAX data. Note that their equa- 
tion (12) agrees with our Eq.(^) within 15 per cent. 

The position of the break in the Lx-SFR relation de- 
pends on the slope of the luminosity function and the value 
of the cut-off luminosity: 



SFRbrcak oc L° 



(13) 



Combined with the slope of the Lx-SFR relation in the low 
SFR regime (Eq.([l2])) and the normalisation of the linear 
dependence in the high SFR limit this opens a possibility to 
constrain the parameters of the luminosity function studying 
the Lx-SFR relation alone, without actually constructing 
luminosity functions, e.g. in distant unresolved galaxies. 



3.5 Lx— SFR relation: comparison with the data 

The solid line in Fig.|^ shows the res ult of the exac t calcu- 



lation of the Lx-SFR relation from Gilfanov et al. (2002) 



The relation was computed for the best fit parameters of the 
HMXB luminosity function determined from the analysis of 
five mostly well studied galaxies from the primary sample 
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Table 4. Sample galaxies from the Hubble Deep Field North and Lynx Field 



Source 


redshift 


^1.4 GHz 


SFR 


So. 5-8 fceV 


Lx 






[^Jy] 


[Mq yr- 1 ] 


[10 -15 erg s _1 cm -2 ] 


[10 40 erg s- 1 ] 


123634.5+621213 


0.458 


233 


28 


0.43 


14.4 


123634.5+621241 


1.219 


230 


213 


0.3 


75.9 


123649.7+621313 


0.475 


49 


8 


0.15 


2.5 


123651.1+621030 


0.410 


95 


9 


0.3 


9.3 


123653.4+621139 


1.275 


66 


69 


0.22 


60.6 


123708.3+621055 


0.423 


45 


4 


0.18 


5.9 


123716.3+621512 


0.232 


187 


5 


0.18 


1.8 


084857.7+445608 


0.622 


320 


71 


1.46 


102 



For two galax ies. 123634.5+621213 and 123651.1+621030, there exist stellar mass estimates obtained with the method of Brinchmann 
& EllisQ pOOOJ) of 4.2 • 10 11 Mq and 7 • 10 10 Mq respectively, which show that the galaxies are dominated by HMXBs (J. Brinchmann, 
private communication). 



(section i.l and Eq. (J6|) ) . Note, that due to the skewness of 
the probability distribution for L to tai in the non-linear, low 
SFR regime the theoretical probability to find a galaxy be- 
low the most probable value (the solid curve in Figj?j) is 
~ 12 — 16 per cent at SFR = 0.2-1.5 Mq yr" 1 and increases 
to « 30 per cent at SFR = 4-5 Mq yr -1 , near the break 
of the Lx-SFR relation. In the linear regime (SFR > 10 
Mq yr -1 ) it asymptotically approaches ~ 50percent, as ex- 
pected. The shaded area around the solid curve corresponds 
to the 68 per cent confidence level including both intrinsic 
variance of the Lx-SFR relation and uncertainty of the best 
fit parameters of the HMXB luminosity function (Eq. (^|) ) . 

Fig.jj] demonstrates sufficiently good agreement between 
the data and the theoretical Lx-SFR relation. Importantly, 
the predicted relation agrees with the data both in the high 
and low SFR regime, thus showing that the data, including 
the high redshift galaxies from Hubble Deep Field North (see 
the following subsection), are consistent with the HMXB 
luminosity function parameters, derived from significantly 
fewer galaxies than plotted in FigQ 

The existence of the linear part at SFR > 5-10 Mq yr -1 
gives an independent confirmation of the reality o f th e cut- 
off in the luminosity function of HMXBs (cf. Sec. 3.2). The 



2001). We cross-correlated the list of the X-ray sources de- 



position of the break and normalisation of the linear part in 
the Lx-SFR relation confirms that the maximum luminos- 
ity of the HMXB sources (cut-off in the HMXB luminosity 
function) is of the o rder of L c ~ 10 40 — 10 41 erg s _1 (see 



Gilfanov et al 



(2002) for more details). Despite the number 
of theoretical ideas being discussed, the exact reason for the 
cut-off in the HMXB luminosity function is not clear and 
significant variations of L c among galaxies, related or not 
to the galactic parameters, such as metalicity or star forma- 
tion rate can not be excluded a priori. However, significant 
variations in L c from galaxy to galaxy would result in large 
dispersion in the break position and in the linear part of the 
Lx-SFR relation. As such large dispersion is not observed, 
one might conclude that there is no large variation of the 
cut-off luminosity between galaxies and, in particular, there 
is no strong dependence of the cut-off luminosity on SFR. 



3.6 Hubble Deep Field North 

In order to check whether the correlation, which is clearly 
seen from Fig. (?) for nearby galaxies, holds for more dis- 
tant galaxies as well we used the data of th e CHANDRA 
observation of the Hubble Deep Field North (Brandt et al 



tected by CHANDRA with t he catalogue o f radio sources 



detected by VLA at 1.4 GHz ([Richards 



identifications of Richards et al 



2000). Using optical 



(1998) and redshifts from 



Cohen et al. (200C) we compiled a list of galaxies detected 



by CHAND RA and classifi e d as s piral or irregular /merger 
galaxies by Richards et al. (1998) and not known to show 
AGN activity. The K-correction for radio luminosity was 
done assuming a po wer law sp e ctrum and using the radio 
spectral indices from Richards (2000). The X-ray luminos- 



ity was K-corrected and transfor med to the 2-10 keV e nergy 
range using photon indices from Brandt et al. (J2001). The 
final list of galaxies selected is given in Table An addi- 
tional data point, X-ray flux and r edshift, is ta k en fro m the 
observation of the Lynx Field by stern et al. (2002) 



The 

radio flux is o btained from a cross-correlation of the X-ray 
positions with Oort (1987). 

The star formation rates were calculated assuming that 
the non-thermal synchrotron emission due to electrons accel- 
erated in supernovae dominates the observe d 1.4 GHz lum i- 
nosity and using the following relation from Condon] (1992): 



SFRr 



= 1.9 • 10~ 29 ■ v% Hz ■ L„[erg s" 1 Hz" 



(14) 



where a is the slope of the non-thermal radio emission. 

The galaxies from HDF North and Lynx are shown in 
Fig.[?j by open circles. A sufficiently good agreement with the 
theoretical Lx-SFR relation is obvious. 



4 DISCUSSION 

4.1 Neutron stars, stellar mass black holes and 
intermediate mass black holes 

Two well known and one possible types of accreting objects 
should contribute to the X-ray luminosity function of sources 
in star forming galaxies: 

(i) neutron stars (M ~ 1.4 M ), 

(ii) stellar mass black holes (3 < M/Mq < 20) born due 
to collapse of high mass stars, and 

(iii) intermediate mass (50 < M/M e < 10 5 ) black holes of 
unknown origin. 

Each class of accreting objects is expected to have a max- 
imum possible luminosity, close or exceeding by a factor of 
several the corresponding Eddington luminosity. In a gen- 
eral case we should expect that each of these three types of 
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accreting objects should have its own luminosity function de- 
pending on the mass distribution inside each class (more nar- 
row for NSs, more broad for BHs and probably very broad 
for intermediate mass BHs), properties of the binary and 
mass loss type and rate from the normal star. Therefore, 
the combined luminosity function of a galaxy, containing all 
three types of objects should have several breaks or steps 
(see Fig. ^) which are not present in Fig. ^. Such breaks 
should be connected with the fact that, for example, below 
the Eddington limit for a NS (or at somewhat higher lumi- 
nosity) more abundant NS X-ray binaries might dominate in 
the number of objects, whereas at higher luminosities only 
black holes should contribute due to their higher masses and 
broader mass distribution. Until now CHANDRA data did 
not show any evidence for a break in the luminosity function 
expected in the vicinity or above of Eddington luminosity 
for NS mass. However, such a break must exist, the only 
question is how pronounced and broad it is. 

It is believed that stars with masses higher than 60-100 
Mq are unstable. Therefore there should be an upper limit 
on the mass of BHs born as a result of stellar collapse. Until 
now the most massive known stellar ma ss BH in our Galaxy , 
GRS 1915+105, has amass of ~15 M Q ( |Greiner et ai]^00l| ). 
It is natural that the Eddington luminosity of these objects, 
amplified several times by angular distribution of radiation 
and chemical abundance effects, should result in the maxi- 
mum luminosity of X-ray sources of this type. It is impor- 
tant to mention that 3 years of RXTE/ASM observations 
revealed from time to time super-Eddington luminosities of 
s ome Galactic X-ra y binaries on the level of 3-12 LEdd,NS 
( |Grimm et ai] |2002j ). 

The hypothetical intermediate mass BHs, probably 
reaching masses of ~ 1O 2_5 M , might be associated with 
extremely high star formation rates (BHs merging in dense 
stellar cluster etc.) and are expected to be significantly less 
frequent than ~stellar mass BHs. Therefore the transition 
from the ~stellar mass BH HMXB luminosity function to 
intermediate mass BHs should be visible in the cumulative 
luminosity function. Merging BHs are one possible way of 
rapid growth of super-massive BHs that exist in practically 
all galaxies. To accrete efficiently intermediate mass BHs 
should form close binary systems with normal stars or be in 
dense molecular clouds. 

If the cut-off in the luminosity function, observed at ~ 
few 10 40 erg s^ 1 corresponds to the maximum possible lumi- 
nosity of ~stellar mass BHs and if at L > L c the population 
of hypothetical intermediate mass BHs emerges, it should 
lead to a drastic change in the slope of the _Lx_-SFR relation 
at extreme values of SFR (Gilfanov et al. 2002). Therefore, 
observations of distant star forming galaxies with very high 



Luminosity | erg/ s 
10 37 10 3B 10 39 1 40 



SFR m ght be one of the best and easiest ways to probe the 
population ot intermediate mass black notes. 

4-1.1 Three component luminosity function 

In Fig. |^ we present the result of a simple picture of what 
type of universal luminosity function a very simple model 
of HMXB population synthesis could produce. This picture 
is obvi ously oversimplified but we present it here to show 




- 0.01 
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Figure 8. The upper main figure shows the contributions of neu- 
tron stars (thin solid line), stellar (dashed line) and intermediate 
(dot-dashed line) mass BHs to the differential luminosity func- 
tion. The thick grey solid line is the total differential luminos- 
ity function. For details of the parameters see discussion in the 
text. The figure in the lower left shows the luminosity around 
the Eddington limit for a NS. The luminosity functions shown 
include the simplest assumption that all systems with accretion 
rates above Eddington radiate at the Eddington luminosity (bot- 
tom), two effects allowing super-Eddington luminosities (middle), 
and additionally a 20 per cent uncertainty in the distance esti- 
mate (top). The curves are scaled for clarification. The figure in 
the lower right shows the luminosity around the Eddington lim- 
its for 3-20 Mq BHs. The luminosity functions shown include no 
effect (bottom), and two effects allowing super-Eddington lumi- 
nosities (middle) and additionally a 20 per cent uncertainty in the 
distance estimate (top). The dotted lines denote the uncertainty 
due to SFR of a factor of 2. 



The initial set-up includes parameterisation of the mass 
distributions of NSs and BHs, the distribution of mass trans- 
fer rates in binary systems, and a prescription for the con- 
version of mass transfer rates to X-ray luminosities. 

The probability distribution of NS masses was chosen 
to be a Gaussian distribution with a mean of 1.4 Mq and a 
a of 0.2 Mq. The mass distribution of BHs was chosen to be 
a power law with a slope of 1.1. These numbers ar e similar 
to results of the oretical computations performed by Fryer & 



Kalogera (2001). The mass distribution for BHs is bimodal, 
for stellar mass black holes it ranges from 3-20 Mq, and 
secondly, we include intermediate mass BHs ranging from 
10 2 — 10 5 Mq. We made the simple assumption that their 
mass distribution has the same slope as for stellar mass BHs. 

Normalisations for the probability distributions were 
chosen such that the number of stellar mass BHs is a factor 
of 20 smaller than the number of NSs. This is roughly the 
ratio observed for HMXBs in our Galaxy ( Portegies Zwart & 



that th? simple picture cannot reproduce the smooth lumi- 



Yungclson 



1995 



[ben et al 



nosity function we get from CHANDRA observations of star 
forming galaxies. 



|l995| ; |Grimm et al.| ^002[ ). How- 
ever the ratio of stars with M > 25Mq, BH progenitors, to 
stars with 25M > M > 8M , NS progenitors, is close to 
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1/2 according to the Salpeter IMF. Therefore in principle 
the stellar mass BH curve in Fig. ^ might be much closer 
to the NS curve. The number of intermediate mass BHs is 
assumed to be a factor of 100 less than the number of stellar 
mass BHs in HMXBs. 

The probability distribution of mass transfer rates in 
binary systems is set to be a power law with a slope of -1.6, 
reproducing the observed luminosity function of HMXBs as- 
suming a linear relation between luminosity and mass ac- 
cretion rate. The limits are 0.1 to 10 7 in units of 10 16 g 
s _1 . Mass transfer was assumed to be conservative over the 
whole range, i.e. no mass is lost from the system except for 
super-Eddington sources and wind accretion. The formulae 
for conversion of mass accretion rate to X-ray luminosity are 



ccretion ' C 



(15) 



where n = 0.1 for BHs and r\ = 0.15 for NSs. The mass 
loss rate from the normal star has no strict limit, however 
the X-ray luminosity reaches a maximum at the Eddington 
luminosity and objects with much higher mass accretion rate 
will end up at the Eddington luminosity introducing a peak 
in the luminosity function. 

For illustration we present two sub-figures in Fig. ^| to 
show the evolution from sharp features to a smoother curve 
with the introduction of smearing effects on the luminos- 
ity which is shown in the main part of the figure. The first 
effect is He-accretion when the HMXB is fed by a helium 
rich star which we take to be the case in about 10 per 
cent of the sources. Secondly, in the case of BHs a quasi- 
f lat accretion disk with an electron scattering atmosphere 
(sobolev 1949j ; Chandrasckhar 1950) radiates according to 
(1 + 2.08 ■ cos(i)) cos(i) where i is the inclination angle, pro- 
ducing 2.6 times higher flux i n the direction perpendic ul ar to 



Titarchuk 



the dis k plane than avera ge (3hakura& Sunyaev 1973). 3un 
yaev 



(1985) confirmed that this ratio is similar 
or higher for r adiation comptonized in t he accretion disk. 
For slim disks (Paczynsky & Wiita 1980| ) this ratio should 
be even higher. Moreover to demonstrate the influence of 
distance uncertainties we assumed a variation in distances 
of 20 per cent. All these effects together give a considerably 
smoother curve and permit up to 6 times higher luminosi- 
ties. 

These are only the most simple effects that permit to 
surpass the Eddington limit. Of c ourse other more sop his- 
ticated models lik e jet emissi on fording et al.| |2002[ ) or 
beamed emission ( King ct al.| 2001) or models ta king into 
account stron g magnetic fields in X-ray pulsars (Basko & 
Sunyaefl |l976l ) also can be employed to explain the observed 



luminosity function. 



4-1-2 Wind driven accreting systems 

Our experience with HMXBs in our Galaxy and LMC shows 
that in many sources accretion happens via capture from a 
strong stellar wind (Cen X-3, Cyg X-l, 4U 1700+37, 4U 
0900-40, and possibly SMC X-l, LMC X-l and LMC X- 
4) As we see the majority of Galactic HMXBs are fed by 
stellar wind accretion. There is a very important difference 
between wind accretion onto NSs and BHs. The capture 
radius, r cap t ure = 2G ^ , is proportional to the mass of the 
accreting object and therefore in similar systems BHs should 



have M 2 times larger accretion rates than NSs for the same 
wind parameters. The dependence of the Roche geometry on 
the mass ratio make the dependence on Mbh a little weaker. 



Ivlcapture OC IVlwind ' \ „ , ) , 

Mns 



(16) 



where (3 is between 1.5 and 2. This reason might increase 
the relative BH contribution to the luminosity function in 
star forming galaxies. It is important that 

M ca pture ■ Tj ■ C 2 ■ ,,/5-l ^ 

oc M wlnd ■ Mg H . (17) 

For /3 > 1 it is preferable for BHs to have higher luminosities 
than for NSs. 



4-1.3 Comparison of simulated and observed luminosity 
function 

The discrepancy between the observed luminosity function 
in the right panel of Fig. ^ and our simple model in Fig. ^ is 
obvious. We do not see features in the observed differential 
luminosity function in the vicinity of L^dd for NSs, neither 
a peak &k ~ 0(2) nor a sharp decline at L > Lsdd as 
in the model luminosity function. Furthermore our model 
luminosity function lacks sources in the luminosity range 
10 39 — 10 40 erg s _1 . It seems we should assume that accreting 
stellar mass BHs in star forming regions are more abundant 
than in the Milky Way. 

It is important to note that having all our corrections we 
are getting objects close to the limit of maximum luminosity 
of the observed luminosity functions. 

In Fig. [s] is plotted the total accretion luminosity 
whereas CHANDRA observes only in the range from 1-10 
keV. However X-ray pulsars emit the bulk of their lumi- 
nosity in the range from 20-40 keV. This effect may further 
decrease the importance of the peak at 2- 10 38 erg s _1 . Since 
in elliptical galaxies old X-ray binaries with weak magnetic 
fields, thus having much softer spectra than X-ray pulsars, 
should dominate the population one should expect the im- 
portance of the peak to be larger in ellipticals. 

Our simple analysis demonstrates how difficult it is to 
construct a very smooth luminosity function with the same 
slope over a broad luminosity range, 10 35 -10 40 erg s"\ and 
without sharp features in the vicinity of Eddington luminosi- 
ties. Because so many different processes are involved in dif- 
ferent parts of this huge luminosity range. Our universal lu- 
minosity function based on CHANDRA, ASCA and RXTE 
data has no strong features. The absence of features around 
the Eddington luminosity for NSs should be explained but it 
is also necessary to explain the absence of the abrupt change 
in the luminosity function at higher luminosities when less 
numerous BHs dominate the luminosity function. 

The most obvious shortcomings of this naive model are 
the mass distributions of BHs and NSs, the normalisations 
for BHs, especially for intermediate mass BHs, and the as- 
sumptions of conservative mass transfer and that all super- 
Eddington sources radiate at Eddington luminosity in X- 
rays. It is also very difficult to assume that intermediate 
mass BHs form a continuous mass function with stellar mass 
BHs without a strong break around 20-50 Mq . They should 
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have their own luminosity function with different normal- 
isation and slope. Another problem is connected with the 
formation of binaries with normal stars feeding intermedi- 
ate mass BHs and making them bright X-ray sources. The 
observation of HMXBs in other galaxies will allow to put 
constraints on the combination of these parameters. 

The main concern with with the existence of a feature- 
less universal luminosity function (ULF) is connected with 
the interpretation of the following experimental facts: 

• RXTE/ASM, ASCA and CHANDRA give us informa- 
tion about the low luminosity part of the ULF (Lx < 10 38 
erg s" 1 ) based on the Milky Way, SMC and NGC 1569. 

• CHANDRA data on the other galaxies in Table [l] give 
information about the high luminosity part of the ULF 
(L X >10 3S erg s" 1 ). 

• UV, FIR and radio methods of SFR determination in 
both local and more distant samples of galaxies have signif- 
icant systematic uncertainties, see Table []. 

To resolve these uncertainties arising very close to the Ed- 
dington luminosity for a NS we need to additional data 
permitting to get the slope of the luminosity function in 
Antennae-type galaxies at luminosities significantly below 
10 38 erg s _1 . Furthermore we need to increase the sample of 
nearby galaxies where we can extend the luminosity function 
well above 10 38 erg s _1 . Only this will give full confidence 
that there is no change in the normalisation in the ULF near 
10 38 erg s _1 . 



4.2 Further astrophysically important information 

The good correlation between SFR and total X-ray lumi- 
nosity due to HMXBs and the total number of HMXBs can 
obviously become a powerful and independent way to mea- 
sure SFR in distant galaxies. In addition, this correlation is 
providing us with further astrophysically important infor- 
mation: 

• These data are showing that NSs and BHs are produced 
in star forming regions very efficiently and in very short 
time, confirming the main predictions of stellar evolution. 

• The luminosity function of HMXBs does not seem to 
depend strongly on the trigger of the star formation event 
which might be completely different for the Milky Way and 
e.g. the Antennae where it is the result of tidal interaction 
of two galaxies. 

• The good agreement of the X-ray luminosity - SFR rela- 
tion of HDF galaxies with the theoretical prediction proves 
that the HMXB formation scenario at high redshifts does 
not differ significantly from nearby HMXB formation. 

• The luminosity function provides information that neu- 
tron stars and BHs have a similar distribution of accretion 
rates in all galaxies of the sample available for study today. 

• The luminosity function of HMXBs does not seem to 
depend strongly on the chemical abundances in the host 
galaxy. 

• The existence of well separated X-ray sources is a way 
to look for small satellites of massive galaxies, like SMC. 

The integral X-ray luminosity and X-ray source counts 
are unique sources of information on binaries in distant 
galaxies. Other methods of investigation of SFR (UV, IR, 



radio) rely on the luminosity distribution and number of 
the brightest stars, without a significant dependence on 
the amount of binaries in a high mass star population. On 
the other hand the existence of an observed population of 
HMXBs in another galaxy is possible only in the case if there 
are conditions for formation of close binaries with certain 
mass loss from a normal companion and efficient capture of 
out-flowing stellar wind or Roche lobe overflow by an accret- 
ing object. Detailed observations of X-ray sources in our own 
Galaxy have shown how small the allowed parameter space 
is - this is the reason why the number of X-ray s ources in 
the Galaxy is so small (Illarionov & Sunyaev [l975| ) m com- 
parison with the total number of NSs and BHs and the total 
number of O and B stars. Therefore: 

• The existence of a universal luminosity function of 
HMXBs proves that the formation of close massive X-ray bi- 
naries and their distribution on mass ratio, separation and 
mass exchange rate is similar in all regions of active star 
formation up to redshifts z~l. 



5 CONCLUSION 

Based on CHANDRA and ASCA observations of nearby 
star forming galaxies and RXTE/ASM, ASCA, and MIR- 
KVANT/TTM data on our Galaxy and the Magellanic 
Clouds we studied the relation between star formation and 
the population of high mass X-ray binaries. Within the ac- 
curacy and completeness of the data available at present, we 
conclude that: 

(1) The data are broadly consistent with the assumption 
that in a wide range of star formation rates the luminosity 
distribution of HMXBs in a galaxy can be approximately 
described by a universal luminosity function, whose normal- 
isation is proportional to the SFR (Fig. [j], [|, Although 
the accuracy of this approximation is yet to be determined 
based on a larger galaxy sample and deeper observations, 
we conclude from the rather limited sample available, that 
it might be of the order of ~50 per cent or better. 

In differential form the universal luminosity function can 
be approximated as a power law with a cut-off at L c ~ 2-10 40 
erg s _1 : 



dN 
dL 3 s 



(3.3t£i) • SFR • Lag 1 ' 61 "' 12 for L < L c 



(18) 



where SFR is measured in units of Mq yr 1 and L38 
L/10 38 erg s . In cumulative form it is correspondingly: 



210" 



N(> L 38 ) = (5.4^7) • SFR ■ (L 3 ~8 ±U 

(19) 

Although more subtle effects can not presently be ex- 
cluded (and are likely to exist), we did not find strong non- 
linear dependences of the HMXB luminosity function on 
SFR. We neither found strong dependences of the HMXB 
luminosity function on other parameters of the host galaxy, 
such as metalicity or star formation trigger. 

(2) Both the number and total luminosity of HMXBs in 
a galaxy are directly related to the star formation rate and 
can be used as an independent SFR indicator. 

(3) The total number of HMXBs is directly proportional 
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to SFR (Fig. 
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SFR [M yr-J = (20) 

(4) The dependence of the total X-ray luminosity of a 
galaxy due to HMXBs on SFR has a break at SFR « 4.5 
M yr _1 for M > 8 M Q . 

At sufficiently high values of star formation rate, 
SFR>4.5M Q yr _1 (L 2 _i kc v>3 10 40 ergs" 1 correspond- 
ingly) the X-ray luminosity of a galaxy due to HMXBs is 
directly proportional to SFR (Fig.[ 



SFRIMoyr- 1 ] = ^ keV , 
1 OJ J 6.7 ■ 10 39 erg s" 1 



(21) 



At lower values of the star formation rate, SFR < 4.5 Mq 
yr _1 (La-io kov <3 ■ 10 40 erg s _1 ), the Lx - SFR relation 
is non- linear: (FigM): 



SFR[M yr _1 ] = 



L 



2-10 kcV 



2.6 ■ 10 39 erg s" 1 



(22) 



The non-linear Lx — SFR dependence in the low SFR 
limit is not related to non-linear SFR dependent effects in 
the population of HMXB sources. It is rather caused by non- 
Gaussianity of the probability distribution of the integrated 
luminosity of a population of discrete sources. We will give 
this a more detailed and rigo rous treatment in a forthcoming 
paper (Gilfanov et al. 2002). 



(5) Based on the data of CHANDRA observations of the 
Hubble Deep Field North we showed, that the relation (J2TJ) 
between the SFR and the X-ray luminosity of a galaxy due to 
HMXBs holds for distant star forming galaxies with redshifts 
as high as z = 1.2 (Fig. M). 

(6) The good agreement of high redshift observations 
with theoretical predictions and the fact that X-ray obser- 
vations exclusively rely on the binary nature of the sources 
is evidence that not only the amount of star formation at 
redshifts up to ~1 can be easily obtained from the above re- 
lations but also that the HMXB formation scenario is very 
similar at least up to this redshift. 

(7) The entire existence of the linear regime in the Lx~ 
SFR relation is a direct consequence of the existence of a 
cut-off in the luminosity function. The position of the break 
in the Lx — SFR relation depends on the cut-off luminosity 
L c in the luminosity function of HMXB as SFRb rC ak <x L" ~ , 
where a is the differential slope of the luminosity function. 
Combined with the slope of the Lx~SFR relation in the 
low SFR regime (Eq.fll3)) this opens a possibility to con- 
strain the parameters of the luminosity function studying 
the Lx-SFR relation alone, without actually constructing 
the luminosity functions, e.g. in distant unresolved galaxies. 

Agreement of the predicted Lx — SFR relation with the 
data both in high and low SFR regime (Fig.Q) gives an inde- 
pendent evidence of the existence of a cut-off in the luminos- 
ity function of HMXBs at L c ~ several x 10 40 erg s~ . It also 
indicates that Lx — SFR data, including the high redshift 
galaxies from Hubble Deep Field North, are consistent with 
the HMXB luminosity function parameters, derived from 
significantly fewer galaxies, than plotted in FigQ 
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